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Abstract 

HIV-1 replication is dependent on binding of the viral capsid to the host protein cyclophilin A (CypA). Interference with 
cyclophilin A binding, either by mutations in the HIV-1 capsid protein (CA) or by the drug cyclosporine A (CsA), inhibits HIV- 
1 replication in cell culture. Resistance to CsA is conferred by A92E or G94D substitutions in CA. The mutant viruses are also 
dependent on CsA for their replication. Interestingly, infection of some cell lines by these mutants is enhanced by CsA, while 
infection of others is not affected by the drug. The cells are thus termed nonpermissive and permissive, respectively, for 
infection by CsA-dependent mutants. The mechanistic basis for the cell type dependence is not well understood, but has 
been hypothesized to result from a dominant-acting host factor that blocks HIV-1 infection by a mechanism that requires 
CypA binding to the viral capsid. In an effort to identify a CypA-dependent host restriction factor, we adopted a strategy 
involving comparative gene expression analysis in three permissive and three non-permissive cell types. We ranked the 
genes based on their relative overexpression in non-permissive cell types compared to the permissive cell types. Based on 
specific selection criteria, 26 candidate genes were selected and targeted using siRNA in nonpermissive (HeLa) cells. 
Depletion of none of the selected candidate genes led to the reversal of CsA-dependent phenotype of the A92E mutant. 
Our data suggest that none of the 26 genes tested is responsible for the dependence of the A92E mutant on CsA. Our study 
provides gene expression data that may be useful for future efforts to identify the putative CypA-dependent HIV-1 
restriction factor and in studies of other cell-specific phenotypes. 
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Introduction 

Cyclophilin A (CypA) is a cellular peptidyl-prolyl isomerase that 
interacts with the HIV-1 capsid and is important for productive 
infection by the virus [1-5]. Interaction of CypA with the viral 
Gag polyprotein in the producer cells leads to incorporation of 
CypA in the budding virions [5] , yet it is the interaction of CypA 
with the incoming viral capsid in the target cell that appears to 
account for the role of CypA in HIV- 1 infection [6,7]. CypA binds 
to an exposed loop on the surface of the CA protein [2]. The 
CypA-binding loop consists of Pro85 to Pro93, with Gly89 and 
Pro90 constituting the binding site for CypA [2]. Preventing the 
CA-CypA interaction by the immunosuppressive drug cyclospor- 
ine A, which targets all cyclophilins [8], or by mutating CypA- 
binding residues in CA leads to impaired infectivity in most cell 
types [1,5-7,9-14]. The effect of CypA appears to occur during a 
post-entry step of the virus life cycle subsequent to reverse 
transcription [12-14]. However, the exact mechanism by which 
CypA promotes HIV-1 replication remains unknown. 

Inhibition of the CypA-C A interaction in HeLa cells by CsA, or 
its analogs leads to impaired HIV-1 replication [15-17]. However, 
these inhibitors have minimal effect on the early post- entry steps 



of the virus life cycle in HeLa cells [6,7,18,19]. Passage of HIV-1 
in HeLa-CD4 + cells in the presence of CsA led to selection of CsA- 
resistant CA mutants A92E and G94D [15,20]. These substitu- 
tions do not detectably alter CypA-CA binding interactions [20]. 
Additional mutants exhibiting CsA resistance have also been 
identified outside the CypA binding loop [21-23]. Interestingly, in 
some cell types replication of the CsA-resistant mutants requires 
CsA or depletion of CypA, indicating that CypA inhibits infection 
in these cells [6,7,24]. The mechanistic basis for this cell-specific 
restriction is not well understood. The cell-type specific nature of 
the CsA-dependent mutants suggests the presence of an inhibitory 
factor in non-permissive cells. Alternatively, a host factor might 
facilitate infection of the permissive cell lines by the mutants. To 
distinguish between these two possibilities, Song and Aiken [18] 
generated heterokaryons by fusion of permissive 293T cells and 
nonpermissive HeLa-P4 cells. Infection of the heterokaryons by 
the mutants was enhanced by addition of CsA, suggesting the 
presence of a dominant-acting, CypA-dependent restriction factor 
in cell lines not permissive to infection by the mutants. 

A common characteristic of all the known CsA-dependent 
mutants is that they are rescued by a specific second-site 
suppressor mutation [13,22,23] suggesting marked phenotypic 
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similarity, thus they can be grouped together as a single class. In 
the present study, we sought to identify a cell-specific host factor 
that restricts this class of mutants in a CypA-dependent manner, 
using A92E as a representative example. We quantified expression 
of genes in permissive and non-permissive cell lines using a human 
gene microarray and ranked them in order of fold-expression in 
non-permissive vs. permissive cell lines. We selected 26 genes with 
at least 3-fold higher expression in non-permissive cell lines, 
knocked them down individually in a non-permissive cell line, and 
analyzed the effects on CsA-dependent infection by HIV-1-A92E. 
We did not identify a gene responsible for reduced infectivity of 
A92E mutant in non-permissive cell lines. Nonetheless, the gene 
expression data reported here can be employed to guide future 
efforts to identify the putative CypA-dependent restriction factor 
as well as other host genes responsible for other cell-specific 
phenotypes. 

Methods and Materials 

Plasmids and Chemicals 

A92E CA mutant was subcloned from HIV-1 proviral DNA 
construct R9 [25] by transfer of Apal-Spel or BssHII-Spel 
restriction fragments into HIV-GFP [26], an envelope-defective 
pNL4-3-based HIV-1 reporter virus clone encoding green 
fluorescent protein (GFP) in place of Nef. The presence of the 
mutations in the final constructs was confirmed by DNA 
sequencing. Plasmid pHCMV-G encodes a vesicular stomatitis 
virus G (VSV-G) protein under the control of the human 
cytomegalovirus (CMV) promoter [27]. Cyclosporine (CsA) was 
purchased from Calbiochem and dissolved in ethanol to produce a 
1 mil stock concentration. 

Cells and Viruses 

Jurkat, CEM, H9 and HOS cells were obtained from NIH 
AIDS Research and Reference Reagent Program as contributed 
by Dr. Arthur Weiss, Dr. J.P.Jacobs, Dr. Robert Gallo and Dr. 
Nathaniel Landau, respectively [28-32]. HeLa and 293T cells 
were kind gifts by Dr. Eric Freed and Dr. Inder Verma 
respectively [33,34]. HeLa, HOS and 293T cells were cultured 
in Dulbecco's modified Eagle's medium (DMEM) (Cellgro) 
supplemented with 10% fetal bovine serum (FBS), penicillin 
(50 IU/ml), and streptomycin (50 ug/ml) at 37°C with 5% C0 2 . 
Jurkat, CEM and H9 cells were cultured in RPMI media 
supplemented with 10% fetal bovine serum (FBS), penicillin 
(50 IU/ml), and streptomycin (50 Ug/ml) at 37°C with 5% C0 2 . 
Virus stocks were produced by calcium phosphate transfection of 
293T cells [35]. VSV-G-pseudotyped reporter virus particles were 
produced by cotransfection of 1 5 (ig of wild-type (WT) or mutant 
HIV-GFP plasmid and 5 ug of pHCMV-G plasmid DNA. Two 
days after transfection, culture supernatants were harvested, 
clarified by filtration through 0.45-um-pore-size filters, and frozen 
into aliquots at — 80°C. The CA content of virus stocks was 
quantified by a p24-specific enzyme-linked immunosorbent assay 
(ELISA), as previously described [36]. 

RNA Isolation for Microarray 

Total RNA was isolated from -lxlO 6 HeLa-P4, HOS and 
293T cells, and ~5 x 10 6 Jurkat, CEM and H9 cells seeded the day 
prior. RNA was purified with TRIZOL per the manufacturer's 
protocol. RNA was further purified using RNeasy kit (Qiagen) 
followed by on-column DNasel treatment. The RNA concentra- 
tion was determined by spectrophotometry, and RNA integrity 
was analyzed by agarose gel electrophoresis. 



Genome-Wide Expression Analysis with the Human Gene 
1 .0 ST Array 

In order to target each "gene" for the desired level of probe 
coverage on the Gene 1.0 ST Arrays the boundaries of each gene 
were defined by calculating "gene bounds". Conceptually, gene 
bounds are the projection of all exons for a given gene onto the 
genome. The gene bounds were calculated in a hierarchical 
manner, using better-annotated evidence first (i.e., RefSeq) and 
more speculative content later (i.e., Ensembl predictions). For the 
Human Gene 1.0 ST Array, only RefSeq, Ensembl and putative 
complete CDS mRNA from GenBank were used to generate gene- 
bound annotations. To ensure robust gene-level estimates of 
expressed RNA, a target of 25 probes per gene was selected for 
each gene bound. Preference was given to probes already on the 
Human Exon 1.0 ST Array. In short, probes were selected 
uniformly over the gene structure. Some of the gene bounds were 
covered by less than 25 probes from the Human Exon 1.0 ST 
Array, so additional probes were selected from the genome 
sequence for the gene bounds. In these cases new probes were 
picked uniformly over the gene structure. Around 80 percent of 
the probes on the Human Gene 1 .0 ST Array are also present on 
the Human Exon 1.0 ST Array. To prepare DNA for hybridiza- 
tion to the array, total cellular RNA was random-primed to 
generate double-stranded DNA using Affymetrix WT cDNA 
synthesis kit. After a clean-up step, double-stranded cDNA was 
fragmented and end-labeled with terminal deoxynucleotidyl 
transferase before hybridization overnight onto Vanderbilt Micro- 
array Shared Resources (VMSR) Human Gene LOST array. The 
single array contained 28,869 gene-level probe sets composed of 
764,885 distinct probes, with a mean of 28 probes per gene and a 
median of 26 probes per gene. Various control probe sets such as 
hybridization control probe sets for BioB, BioC, BioD and CreX 
pre-labeled spikes, assay control probe sets for the Dap, Phe, Lys 
and Thr polyA unlabeled spikes, putative exon and intron control 
probe sets from putative constitutive genes for use as pseudo 
positive/negative controls and generic background probes were 
also included in the array. All data were RMA normalized using 
Partek. 

siRNA Knockdown and viral Infectivity Assay 

HeLa-P4 cells were transfected with either SMARTpool 
siRNAs (Dharmacon) or scramble siRNA using siPORT NeoFX 
transfection reagent (Invitrogen) per the manufacturer's protocol. 
As a control for functional knockdown, we transfected a parallel 
culture with a well-characterized siRNA against TNP03, a 
protein on which HIV-1 infection is dependent. 72 hours after 
transfection, cells were trypsinized and replated in 24-well plates at 
a density of 25,000 cells/well. The next day, monolayers were 
inoculated with VSV-G-pseudotyped HIV-GFP-A92E in the 
presence or absence of CsA (5 u.M). Twenty-four hours later, 
the cultures were supplemented with 1 ml of additional medium 
and incubated for another 24 h. Cells were then detached with 
trypsin, fixed in 2% paraformaldehyde, and analyzed for GFP 
expression with an Accuri C6 flow cytometer. 

Results 

Cell-type-dependent Phenotype of the A92E Mutant 

Previously, several studies have demonstrated CsA-dependence 
of the A92E mutant in HeLa and H9 cells [6,24] . In order to 
identify other cell types in which infection by the A92E mutant is 
dependent on CsA, we tested several cell lines for VSV-G- 
pseudotyped-A92E infectivity in the presence or absence of CsA. 
As shown previously [6,24], we observed that CsA enhanced 
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infection of HeLa and H9 cells 6.5-fold and 2.0-fold, respectively 
(Table 1). In addition to H9 cells, we employed another T cell line, 
CEM, in which CsA enhanced A92E infection by 1.8-fold 
(Table 1). In contrast to the non-permissive cell lines, we observed 
that CsA exhibited minimal effects on infection by the A92E 
mutant in HOS and 293T cells. Moreover, in Jurkat cells, CsA 
inhibited infection by approximately 50%. 

Human Gene Array 

To identify genes that may play a role in the CypA-dependent 
restriction of the A92E mutant, we performed Human Gene Array 
analysis of the six cell lines in triplicate. In all, expression of 28,869 
genes was analyzed across these six cell types (complete data set 
shown in Table SI). The cell lines were categorized as 
nonpermissive (HeLa, H9, CEM) and permissive (Jurkat, HOS, 
293T). For each gene, the mean expression in the permissive cell 
group was compared to the corresponding value for the non- 
permissive group. In this way, the fold-change in expression in 
non-permissive cells to permissive cells was determined for each of 
the 28,869 human genes. The genes were then sorted by the ratio 
of the expression in non-permissive cells versus permissive cells 
(Table S2). The top 100 candidate genes with at least 3-fold higher 
expression in non-permissive cells compared to permissive cells 
were selected. This list was further scrutinized based on two 
criteria: (1) Expression of a gene in HeLa cells should be higher 
than in any of the permissive cell types because A92E demon- 
strated the most pronounced CsA-dependent phenotype in HeLa 
cells; and (2) Expression of a gene in any of the permissive cell 
types should be lower than in at least two of the non-permissive 
cell types. Using these two criteria, 26 genes (Table 2) were 
selected for analysis of their effects on HIV-1-A92E infection. 
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Analysis of the Candidate Genes 

To test if expression of any of the selected candidate genes is 
necessary for the observed CsA-dependent phenotype of the A92E 
mutant, each of these genes was individually knocked down in 
HeLa cells using SMARTpool siRNAs. For the assays, HeLa cells 
were chosen because the A92E mutant demonstrated the strongest 
enhancement by CsA in this cell type. A scrambled siRNA 
sequence was used as a control for non-specific effects of siRNA 
transfection. The analysis was performed in several experiments 
with subsets of the siRNAs, each of which also contained the 
scrambled control. Seventy two hours after siRNA transfection, 
cells were harvested and plated for infection. Cells were 
transduced with VSV-G-pseudotyped HIV-GFP-A92E mutant 
particles in the presence or absence of CsA sixteen hours after 
plating. Transfected cells were analyzed for single-cycle infection 
by the mutant forty-eight hours after inoculation. No apparent 
cytotoxicity was observed with knockdown of these candidate 
genes using 200 nM of their respective siRNAs. We reasoned that 
depletion of the restriction factor would result in enhancement of 
HIV-1-A92E infection in the absence of CsA, and would have 
minimal effects on infection in the presence of CsA, resulting in a 
net reduction in the extent to which CsA enhances infection of the 
virus. We observed that depletion of any of the 26 candidate genes 
did not result in complete rescue from CsA dependence (Table 3). 
However, depletion of nine genes including TLL1, AGMAT, 
JVEDD9, PRSS21, IFITM1, PLEKHG1, SLC16A6, KAMI and 
ACSL5 resulted in 1 .6-2.7-fold enhancement of A92E infectivity in 
the absence of CsA and 1.4— 2.0-fold decrease in CsA enhance- 
ment of A92E infection. In addition, depletion of these 26 genes 
had no significant effect on infectivity of the wild-type virus (data 
not shown). 
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Discussion 



The mechanism by which CypA restricts infection by some 
HIV-1 CA mutants is unknown. In this study we sought to test the 
hypothesis that a dominant-acting CypA-dependent restriction 
factor is responsible for CsA-dependent enhancement of these CA 
mutants. Assuming that the putative factor is differentially 
expressed in permissive and nonpermissive cells, we performed 
gene expression analysis on three permissive and three non- 
permissive cell lines. Based on our selection criteria, we focused on 
26 candidate genes and individually knocked them down in HeLa 
cells using siRNAs. In order to thoroughly inhibit gene expression 
while minimizing cytotoxicity, we employed SMARTpool siRNAs 
and transfected cells using three different concentrations of 
siRNAs (8nM, 40 nM and 200 nM; Table 3 and data not shown). 
The concentrations of siRNAs used spanned the range of 
manufacturer-recommended concentrations. We expected a 
reversal of A92E phenotype when the putative restriction factor 
is depleted. Our data showed that depletion of at least nine genes 
resulted in 1.7-2.7-fold increase in A92E infectivity in the absence 
of CsA and 1.4-2.0-fold decrease in CsA dependence. Interest- 
ingly, of the nine genes, IFITM1 was demonstrated as a potent 
inhibitor of influenza A H 1 N 1 , dengue virus, West Nile virus and 
HIV-1 [37,38]. Nevertheless, our results do not provide clear 
evidence for involvement of any of the 26 genes analyzed here in 
CsA enhancement of infection by HIV-1-A92E. It remains 
possible that low knock-down efficiency precluded identification 
of the putative gene. Alternatively, CypA-dependence could 
involve multiple gene products, thus single knock-downs may 
not reveal the restriction factor. Although this study focused on 
differences in expression levels, it is possible that single nucleotide 
polymorphisms could explain the phenotypic differences between 
cell types. Notably, previous studies have demonstrated that 
polymorphisms in the regulatory region of PPIA gene, which 
encodes CypA, influence susceptibility to HIV-1 infection or 
disease progression [39-41]. 

The CA mutants A92E and G94D were originally isolated as 
CypA-independent mutants by serial passage of wild-type HIV-1 
in HeLa cells [15,20]. These mutants are also dependent on CsA 
for their replication. Subsequently, it was established that the CsA- 
dependence of these mutants varies among cell lines [6,7,19,24]. 
These mutants replicate efficiently in the presence of CsA in HeLa 
and H9 cells but not injurkat, HOS, TE671 and 293T cells. We 
identified another human T cell line, CEM, in which the A92E 
mutant requires CsA for infection. The cell-type dependence of 
the CsA-enhanced infection by these mutants is not understood. 
Evidence from heterokaryon analysis suggests that a dominant- 
acting CypA-dependent host factor inhibits A92E infection in 
HeLa cells [18]. Based on higher expression of CypA in non- 
permissive H9 cells compared to that in non-permissive Jurkat 
cells, it was suggested that differing CypA levels could explain the 
cell-type-dependent phenotype of the mutants [24]. Moreover, 
permissive TE671 cells overexpressing CypA were rendered non- 
permissive to infection by the mutants [19]. In contrast, 
overexpression of CypA in permissive 293T cells did not lead to 
enhancement of infection by the mutants upon CsA treatment 
[18]. Li et al. reported that the A92E mutant capsid undergoes 
more rapid uncoating in HeLa cells than injurkat cells, suggesting 
that CypA may inhibit this virus by destabilizing the capsid. [12]. 
Paradoxically, we have previously shown that CypA can stabilize 
wild type HIV-1 cores in vitro [42]. Moreover, a recent study 
demonstrated that CypA-CA interaction inhibited nuclear entry of 
A92E and increased the levels of pelletable CA protein in HeLa 
cells. The authors suggested that a dominant-acting restriction 
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factor restricts the A92E mutant by stabilizing the viral capsid in a 
CypA-dependent manner [43] . Precedence for a CypA-dependent 
restriction factor exists in Old World monkey cells. In these cells, 
binding of CypA to the HIV- 1 capsid leads to diminished HIV- 1 
infectivity by sensitizing the capsid to restriction by TRIM50C [44— 
46] . In human cells, CypA interaction with CA seemed to protect 
HIV-1 from a saturable restriction factor [47]. Initial studies in 
human cells suggested that huTRIM5a could restrict HIV-1 
infection in a CypA-dependent manner [11,48]. However, later 
studies demonstrated that endogenous human TRIM5a modesdy 
inhibits HIV-1 and inhibition of CypA-CA interaction does not 
diminish its antiviral activity [46,49,50]. Collectively, the evidence 
is consistent with the presence of an unknown CypA-dependent 
restriction factor in cell types that are nonpermissive to CsA- 
dependent mutants. Interestingly, three recent studies have 
identified the human MxB protein as a CypA-dependent HIV-1 
restriction factor that is induced by interferon [51-53]. 

One strong candidate for the CypA-dependent HIV-1 restric- 
tion factor is the host protein CPSF6. A C-terminally truncated 
CPSF6 (mCPSF6-358) protein, which localizes to the cytoplasm, 
was demonstrated to bind HIV-1 CA and inhibit viral nuclear 
entry [54]. The N74D substitution in CA renders the virus 
resistant to CPSF6-358 by inhibiting binding of CPSF6-358 to the 
capsid. Like wild-type HIV-1, A92E mutant is inhibited by 
CPSF6-358 [55]. Interestingly, the N74D/A92E double mutant is 
resistant to CPSF6-358 restriction, suggesting that N74D mutation 
rescues A92E presumably by preventing CPSF6 binding [55]. 
Moreover, infection by another CsA-dependent mutant T54A is 
rescued by addition of A105T, which is present in the CA binding 
site for CPSF6 and confers resistance to CPSF6-358 [56]. It was 
recently reported that CsA-dependent mutants are partially 
rescued by depletion of endogenous CPSF6 in HeLa cells [57]. 
However, expression levels of CPSF6 in permissive and non- 
permissive cell lines were similar in our study (Table 4), thus 
variations in its expression do not correlate with the cell- 
dependence of infection by the HIV-1-A92E mutant. Therefore, 
it remains unknown whether the cell-type-dependent phenotype of 
CsA-dependent mutants is accounted for by CPSF6. 

The goal of this study was to identify a novel restriction factor in 
non-permissive cell types. We chose three permissive and three 
non-permissive cell types and performed human genome micro- 
array analysis on these cell types. We grouped cell lines based on 
their permissivity to HIV-1-A92E infection and then compared 
the mean levels of gene expression in permissive versus non- 
permissive cell types to generate a ranked list of candidate genes. 
Our approach was mainly driven by the observation that the 
magnitude of CsA-dependence of A92E was ~ 3-fold higher in 
HeLa cells than in other non-permissive cell types. Due to limited 
resources, we analyzed only 26 candidates from the first 100 genes 
whose expression in non-permissive cell types group was 3-fold or 
higher than in permissive cell types group. A limitation of this 
approach is that it ignores differences in CsA dependence within 
each group. An alternative approach involves correlating the fold- 
difference in CsA-dependent phenotype between cell types with 
fold-change in gene expression. To this end, we calculated a 
Pearson correlation coefficient for each gene to correlate the CsA- 
dependence of the mutant with gene expression (Table S3). This 
approach accounts for differences in CsA-dependence within each 
group. Of the 26 genes selected by the first approach, four genes 
(SLC16A6, MCTP, KAMI and TLL1) had Pearson correlation 
coefficients in the range of 0.97-0.98 and were among the top- 
ranked 1 00 genes by the Pearson method. Moreover, none of the 
selected hits contain a CypA domain or associate with either CypA 
or CypA domain containing protein. Expression of genes encoding 



known HIV-1 capsid interaction partners [4,14,54,58-61] were 
also analyzed (Table 4). Some studies have suggested that CypA 
levels dictate the CsA-dependent phenotype of the A92E mutant 
[19,24]. However, the gene expression data showed a negative 
correlation between CYPA expression and CsA- dependence of the 
A92E mutant. Similarly, expression of CPSF6 was also negatively 
correlated with CsA-dependence of the mutant. NUP153, 
JVUP358, MJP98, ERK and CYPB had positive but poor 
correlation with the CsA-dependent phenotype. Collectively, the 
alternative approach suggests that the known capsid-interacting 
proteins are unlikely to determine CsA-dependence of the A92E 
mutant. 

In summary, although our limited study failed to identify the 
putative restriction factor, the approach using DNA microarray 
and the results thereof can be employed by other interested 
investigators for a more detailed analysis. The gene expression 
data may also be useful to investigators interested in other 
phenotypes that vary among the cell lines employed in this study. 

Supporting Information 

Table SI Human gene array. The table shows RMA 
normalized (log2 scale) values, mRNA accession numbers and 
gene symbols for all the genes analyzed by Human Gene 1 .0 ST 
Array. The gene expression analysis was performed on six cell lines 
in triplicate. 
(RAR) 

Table S2 Fold change in expression of genes in non- 
permissive vs permissive cells. Mean expression values of 
each gene in permissive and non-permissive cell types were 
calculated. The mean value was used to calculate fold change in 
expression of each gene in permissive (CA1-9) vs non-permissive 
(CA10-18) cell types. Fold expression values with a (— ) sign 
indicate higher expression in non-permissive cells, values with a (+) 
sign indicate higher expression in permissive cells. Columns also 
include annotation, uncorrected p-value, two p-values corrected 
by multiple testing correction (Bonferroni which is very stringent 
and then STepup also called Benjamini-Hochberg which controls 
the error in a different way but is less stringent). 
(XLSX) 

Table S3 Correlation between CsA-dependent pheno- 
type of the mutant with gene expression. Pearson 
correlation coefficient for each gene was calculated to correlate 
the CsA-dependence of the mutant (calculated in Table 1) with 
gene expression. The list is sorted by the Pearson correlation 
coefficient value. Slope values were obtained by plotting CsA- 
dependent phenotype (X-axis) in all the cell lines against the gene 
expression values (Y-axis). 
(XLSX) 
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